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I.  INTRODUCTION 

By  contrast  with  the  number  of  papers  devoted  to  the  investigation  of  noise 
radiated  by  a region  of  fully-developed  turbulence,  only  a few  have  addressed  the 
sound  field  generated  by  the  laminar- to- turbulent  transition  zone  [1  - 4). 

Ffowcs  Williams  [1]  discusses  the  transition  zone  noise  in  terms  of  the  initial 
formation  of  turbulence  spots.  The  elegant  analysis  of  a three-dimensional 
disturbance  growing  in  accordance  with  the  equations  of  linear  stability  theory 
by  Brooke  Benjamin  [5]  was  incorporated  into  Lighthill's  theory  [6]  for  the 
radiating  component  of  flow  noise.  A quadrupole  source  was  assumed  and  it  was 
thus  shown  that  sound  radiation  would  occur  essentially  only  when  the  wall  was 
of  pressure  releas.e  type.  For  an  assumed  rigid  surface  the  imaging  effect 
transforms  the  quadrupoles  into  much  less  efficient  octupole  radiators.  Further- 


more, even  for  the  soft  surface,  the  analysis  shows  that  it  is  only  at  the  very 
beginning  of  spot  formation  where  sound  is  generated.  Farther  downstream  the  spot 
grows  exponentially,  but  the  sound  pressure  is  found  to  decrease  exponentially  with 
time.  Of  course  it  could  be  argued  that  the  linear  stability  model  breaks  down  as 
a spot  develops  which  may  account  for  this  conclusion.  Ffowcs  Williams  points  out 
that  the  theory  may  be  applicable  only  if  sound  is  truly  generated  in  the  early 
stages  of  instability,  but  late  enough  such  that  the  asymptotic  theory  of  Brooke 
Benjamin  is  established,  while  at  the  same  time  the  spot  amplitude  is  still  small. 

Dolgova  [4]  considered  the  planar  flow  over  a rigid  surface  and  analyzed  the 
sound  radiated  by  only  the  Tollmien  - Schlichting  wave  growth.  As  in  linear 
stability  models  she  let  a pressure  wave  grow  exponentially  in  the  direction  of 
flow,  and  then  applied  Fourier  and  Hankel  transforms  to  obtain  an  expression  for 
the  sound  pressure  and  directivity  function.  At  low  frequencies  the  pattern  is 
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dipole  while  at  higher  frequencies  the  main  lobe  tilts  to  the  downstream  direction, 
gets  sharper,  and  develops  several  side  lobes.  The  resulting  expression  for  the 
sound  pressure,  however,  shows  no  Mach  number  dependence. 

Natural  transition  is  characterized  by  essentially  three  distinct,  flow 
regimes.  In  the  early  stages  of  instability,  the  laminar  boundary  layer  becomes 
disturbed  in  a linear,  wavelike  manner.  It  is  here  where  the  theory  of  Dolgova  [4] 
is  applicable.  Farther  downstream,  these  linear  disturbances  become  more  non- 
linear and  three  dimensional.  The  analysis  of  Ref.  [1]  appears  to  be  applicable 
in  this  regime,  but  breaks  down  rapidly  as  turbulent  bursts  begin  to  form.  Within 
the  third  flow  regime,  where  the  boundary  layer  intermittently  alternates  between 
laminar  and  turbulent,  is  where  no  fundamental  flow  noise  theory  has  yet  been 
developed.  We  hypothesize  that  this  region  of  intermittent  boundary-layer  flow 
may  give  rise  to  sound  radiation  [2,  3].  We  know  from  elementary  boundary-layer 
theory  that  the  mean  velocity  gradients  and  profiles  are  quite  different  between 
laminar  and  turbulent  flow  over  surfaces.  In  intermittent  flow,  we  would  expect 
that  the  wall  shear  stress  and  the  mean  boundary  layer  velocities,  both  parallel 
and  normal  to  the  surface,  undergo  gross  fluctuations  in  time.  According  to 
Lighthill's  analogy  [6],  fluctuating  velocities  give  rise  to  quadrupole  noise 
sources,  while  fluctuating  wall  shear  stresses  give  rise  to  dipole  sources  [7]. 

At  very  low  Mach  numbers,  M,  this  latter  source  would  be  expected  to  dominate  the 

noise  field.  This  is  deduced  from  the  fact  that  the  acoustic  efficiency  of  dipole 

3 5 

sources  is  of  the  order  M , while  that  of  quadrupole  sources  is  of  order  M , viz., 

Ross  [8]. 

In  this  paper  we  will  approach  analytically  the  problem  of  noise  generated  by 
intermittent  boundary-layer  flow  as  it  occurs  in  natural  transition  on  a flat  plate 
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We  will  assume  the  plate  to  be  infinite  in  extent  and  to  be  acoustically  rigid.  We 
will  further  assume  that  incompressible  hydrodynamics  will  describe  the  transitional 
process.  Our  goal  is  to  derive  expressions  for  the  radiated  noise  spectrum  and 
acoustic  efficiency.  The  spectrum,  because  of  our  infinite  extent  assumption,  will 
be  derived  in  terms  of  the  sound  pressure  radiated  per  unit  spanwise  width  of 
transition.  A comparison  of  predicted  spectra  with  experimental  data  obtained  on 
the  surface  of  a buoyantly-propelled  vehicle  is  also  presented. 

II.  ANALYSIS 

Consider  the  uniform  flow  over  a flat  plate.  If  a laminar  boundary  layer 
begins  to  form  at  the  origin  of  our  Cartesian  coordinate  system,  and  if  the  flow 
direction  is  in  the  x -direction,  then  the  flow  will  ultimately  become  nonlinearly 
unstable  at  the  downstream  line  x^=Xq  (see  Fig.  1).  There  is  then  a short  distance. 
Ax,  over  which  turbulent  "bursts"  occur.  Along  x^  = x^  + Ax,  the  flow  becomes 
fully  turbulent  while  for  x^  < Xq,  the  flow  is  assumed  completely  laminar.  In  the 
analysis  to  follow,  we  will  treat  SAxdx^  as  our  source  volume,  where  5 is  the 
boundary  layer  thickness. 

The  equation  that  describes  the  sound  radiation  from  fluid-dynamic  sources  is 
the  well-known  Lighthill  equation  [6]: 


at2  ax  2 3xi9xj 


(i) 


where  p'  is  the  fluctuating  density,  t is 
medium  surrounding  the  source  region,  and 
by: 


time,  c is  the  sound  velocity  in  the 
T^  is  Lighthill' s stress  tensor  given 
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Tij  = Puiuj  - CTlj  + (p,"p,C  )<Sij 


(2) 


where  p is  the  mean  density. 


Composing  this  tensor  is  pu^u^ , the  fluctuating  Reynolds  stress  tensor;  0^ , 


2 

the  fluctuating  viscous  shear  stress  tensor;  and  (p'-p'c  ),  a term  that  relates  to 
heat  conduction  or  nonlinearity.  For  hydrodynamic  flows,  this  last  term  is  zero 
which  means  that  the  fluctuating  density  is  related  to  the  fluctuating  pressure  by 


p’  - P’c2 


(3) 


From  this  point  on,  we  will  drop  the  prime  from  the  fluctuating  quantities,  and 
solve  Eq.  (1)  for  the  acoustic  pressure  through  use  of  Eq.  (3). 


A.  for  Boundary-Layer  Transition 


Within  the  region  x^  £ x^  + Ax,  the  flow  intermittently  changes  from 

laminar  to  turbulent  regimes.  A typical  turbulence  burst  is  shown  schematically  in 
Fig.  2.  A burst  will  grow  as  it  travels  downstream  at  a mean  convection  velocity, 
uc«  The  locus  of  a typical  burst  forms  a wedge  of  apex  angle  2a.  Emmons  [9]  first 
observed  a to  be  on  the  order  of  9.6°  for  flat  plate  flows.  This  was  reconfirmed 
later  by  Schubauer  and  Klebanoff  [10],  although  it  should  be  pointed  out  that 
Farabee,  et.al.  [11]  found  propagation  angles  significantly  greater  than  this. 

In  reference  to  the  elevation  view  of  Fig.  2,  we  might  anticipate  that  at  a 


given  point  on  the  x^-axis,  the  mean  velocity  profile  and  the  wall  shear  stress 


will  undergo  gross  fluctuations  in  time  as  bursts  are  swept  by.  We  will  assume 

that  these  fluctuations  are  much  larger  than  the  velocity  and  shear  stress  fluctuations 

that  occur  within  a burst  itself  or  within  the  fully-developed  turbulent  boundary 


layer  that  occurs  at  x^  > Xq  + Ax.  For  example,  the  wall  shear  stress  under  a 


turbulent  boundary  layer  can  be  calculated  from  [12]: 


J 
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TT^X1^  = °-0288Puo  ^U0X1^  ’ 


while  at  some  other  instant  of  time,  when  the  boundary  layer  is  laminar  at  x^: 


Tl(Xi>  = 0. 332pvuQ(vx1/up)  , 


where  V is  the  kinematic  viscosity.  Thus,  as  the  flow  regime  alternates  between 
laminar  and  turbulent,  we  might  expect  fluctuations  in  the  wall  shear  stress  on  the 
order  of  o(x^),  where 


0(xi>  _ 

tt(xi)  = tt 


-3/10 


1 - 11. 53 (Re  ) 
X1 


and  Re^  = u^x^/v  which  is  the  local  value  of  length  Reynolds  number.  Typically, 

Re  - 3 x 10  , so  cx(xn)  - 0.9tt(x_).  Although  not  quantitatively  proven,  we  would 
x0  u v 

expect  that  fluctuations  of  this  magnitude  far  exceed  the  wall  shear  stress  fluctuations 
that  may  occur  beneath  a fully-developed  turbulent  boundary  layer  (compare,  for 
example.  Fig.  16.17  with  Fig.  18.4  of  Schlichting  [12]). 

Consequently,  Eq.  (2)  for  T may  be  written 


where 


Tij  ■ cVj  ■ W0  • 


W'>  - °12|x  - 0 ~ VljMM 
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which  is  the  fluctuating  wall  shear  stress  due  to  the  creation  and  convection  of 
turbulence  bursts.  If  we  further  assume  (because  of  the  slow  rate  at  which  a burst 


spreads  laterally)  that  the  mean  velocity  fluctuation  in  the  x^-direction  is  small 
compared  to  u^  and  U2,  then  will  have  only  three  dominant  components  in  addition 
to  Tg^.t). 


We  have  not  yet  been  explicit  regarding  the  time  dependence  of  t^,  u^,  and  U£. 
It  will  prove  expedient  first  to  solve  Eq.  (1)  subject  to  our  assumed  form  of 
and  the  appropriate  acoustic  boundary  conditions. 


B.  Formal  Solution 

The  solution  of  Eq.  (1)  for  flow  over  arbitrary  surfaces  was  first  constructed 
by  Curie  [13].  Powell  [14]  showed  that  when  the  surface  is  rigid  and  planar,  one 
could  consider  a new  extended  flow  field  obtained  by  reflection  of  the  original  one 
in  the  plane  X2  = 0,  as  illustrated  in  Fig.  3(a).  Using  Powell's  result,  we  find 


where  the  integrands  are  functions  of  the  retarded  time  variable  (t-r/c)  and  primes 
refer  to  the  image-flow  side  of  the  plane.  The  radial  coordinate  defining  the 
point  of  observation  is  r = |x  - f||>  where  r|  denote  the  source  volume  coordinates. 
By  making  the  far-field  assumption,  we  can  let  r - r'  = |x|,  and  then  by  making  the 
standard  transformations  from  spatial  derivatives  to  time  derivatives,  Eq.  (9) 


becomes: 
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(10) 


Each  of  the  terms  composing  Eq.  (10)  represent  a hydrodynamically  - generated 

acoustic  source.  The  direction  cosine  preceeding  each  integral  tells  us  the  order 

and  orientation  of  each  of  these  sources.  We  see  that  the  first  and  fourth  terms 

2 

have  a cos  0 directivity  characteristic  (see  Fig.  1 for  definition  of  0),  indicating 

that  they  are  longitudinal  quadrupoles  with  axes  parallel  to  the  surface.  These 

two  sources,  one  being  in  the  image  flow  combine  constructively  to  form  a single 

longitudinal  quadrupole  of  twice  the  strength  of  one  alone.  The  second  and  fifth 

terms  represent  lateral  quadrupoles  because  of  the  cross-term  products  and  directivity 

function  sin  0 cos  0.  However,  when  combining  a lateral  quadrupole  with  its  image, 

a less  efficient  octupole  results.  We  can  safely  neglect  those  source  terms  involving 

u^  U£  and  u^'u2'  un^er  our  low  Mach  number  assumption.  The  third  and  sixth  terms 

2 2 

also  represent  longitudinal  quadrupoles,  but  because  = sin  0 their  axes  are 

perpendicular  to  the  surface;  when  combined,  another  inefficient  octupole  results  [15]. 
The  last  term  of  Eq.  (10)  includes  the  image  contribution  and  describes  a dipole 
source  with  axis  parallel  to  the  surface.  Fig.  3(b)  illustrates  in  an  elementary 
manner  the  combination  of  these  various  sources. 
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Our  solution  to  Eq.  (1)  thus  reduces  to  the  sum  of  a volume  integral  describing 
longitudinal  quadrupoles  and  a surface  integral  describing  dipoles,  i.e., 

2 2 f 

, pcos  0 3 

P(r,t)  = = 2 

2TTrc  3t 

V 

• • 

T0(oi,t-r/c)dS(Tli)  . (11) 

s 

As  pointed  out  in  the  INTRODUCTION,  we  can  expect  the  dipole  contribution  of  Eq.  (11) 

-2 

to  be  of  the  order  M more  efficient  than  the  quadrupole  contribution.  As  a first 
approximation,  we  would  like  to  accept  this  supposition,  and  examine  in  detail 


cos6  3 
+ 2irrc  3t 


ui  (ni,t-r/c)dvcni) 


P(r,t) 


cos9  3 
27rrc  3t  | 

S 


T0(ni>t-r/c)ds(ni) 


(12) 


when  M « 1. 

C.  Proposed  Model  for  the  Wall  Shear  Stress  Fluctuations 

We  have  developed  the  solution  given  by  Eq.  (12)  upon  the  notion  that  it  is 
the  intermittency  of  the  boundary  layer  within  the  transition  zone  that  creates 
noise.  Intermittent  boundary-layer  flow  is  a situation  where  we  can  construct  a 
mathematical  model  for  the  fluctuating  physical  parameters  of  interest.  This  model 
is  analogous  to  that  which  the  experimentalist  uses  to  distinguish  time  intervals 
when  his  sensor  is  in  irrotational  fluid  from  those  when  it  is  in  turbulent  fluid. 
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I 
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Clearly,  this  is  a zero-one  function  which  we  shall  call  the  indicator  function, 
I(x^,t).  The  use  of  such  a function  to  describe  intermittently  turbulent  flows  is 
not  new,  viz.,  Refs.  [16-18],  among  others.  However,  as  will  become  apparent  below, 
we  will  require  the  space-time  correlation  functions  of  ICx^.t),  and  these  functions 
have  not  been  so  well  investigated  for  boundary-layer  transition  [19]. 

Typically,  the  indicator  function  at  a given  point,  x^,  in  the  source  volume 
may  be  illustrated  as  in  Fig. 4.  This  function  is  zero  when  the  flow  is  laminar  and 
is  unity  when  turbulent.  Also  shown  in  this  illustration  is  the  first  time  derivative 
of  I(x^,t)  which  is  a random  sequence  of  alternating  Delta  functions.  Because 
Eq.  (12)  needs  only  to  be  evaluated  on  the  surface,  the  functional  dependence  of  I 
on  x^  is  not  required.  It  is  reasonable  to  let  I be  statistically  homogeneous  in 
x^  because  of  our  infinite  plane  surface  assumption;  however,  I must  necessarily  be 
non-homogeneous  in  x^. 

The  non-homogeneity  of  I in  the  x^-direction  is  due  to  the  fact  that  the 
boundary  layer  ultimately  changes  from  fully  laminar  to  fully  turbulent.  As  x^ 
increases  beyond  x^  more  and  more  impulse  functions  fill  in  the  time  scale.  The 
time-average  value  of  I(x^,t)  is  appropriately  called  the  intermittency  factor,  i.e.. 


and  represents  the  fraction  of  time  that  the  flow  is  turbulent.  This  is  a pre- 
dictable function  for  boundary-layer  transition  on  flat  plates  [9]  as  well  as  for 
axisymmetric  bodies  [20].  Another  important  mean  property  of  I(x^,t)  is  the 
"burst  frequency,"  N(x^).  This  function  describes  the  expected  number  of  bursts  that 
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occur  at  a given  point  per  unit  time.  It  too,  may  be  predicted  for  most  hydro- 
dynamic  flows  [9,  11,  20].  Fig.  5 shows,  in  principle,  how  the  indicator  function 
varies  in  time  and  space;  also  shown  are  typical  distributions  of  y and  N. 

With  these  definitions  for  I(x^,t),  and  its  mean  properties,  we  propose  that 

T0^Xi* = Ti/Xi^  + Ti^Xi^  0-4) 

where  i may  take  on  the  values  1 and  3.  By  taking  the  time  derivative  of  Eq.  (12) 
inside  the  integrals,  we  find: 

TT  = <Tt'TL)  S = °<Xi)  i(xi>,:)  ‘ <15) 

where  0 was  first  defined  in  Eq.  (6). 

As  this  point  it  is  important  to  emphasize  that  Eq.  (14)  implies  that  the  wall 
shear  stress  is  capable  of  changing  from  a laminar  value  to  a turbulent  value  in  an 
infinitesimally-short  period  of  time  (I's  are  Delta  functions).  Obviously,  this 
would  seem  physically  impossible.  However,  we  do  know,  from  oscillograph  traces  of 
the  velocity  and/or  pressure  fluctuations  that  occur  in  transition  flows,  that  the 
flow  state  can  change  in  extremely  short  periods  of  time  as  turbulent  bursts  are 
created  at  or  swept  by  the  measuring  sensor  [10-12].  Eq.  (14),  although  idealistic, 
may  be  quite  adequate  for  estimating  the  low-frequency  portion  of  the  radiated 
noise  spectrum.  The  techniques  developed  by  Schottky  [21]  in  his  classic  analysis 
of  shot  noise  would  seem  appropriate  here. 

As  a final  point  regarding  Eq.  (14),  we  note  that  its  time  average  is  simply: 

- [l-y(x1)]  tL(x1)  + y(xx)  t^x^ 


(16) 
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which  is  identical  to  that  suggested  by  Emmons  [9]  for  the  transition  region. 
Several  experiments  have  verified  this  equation,  e.g.,  Dhawan  and  Narasimha  [22]. 

D.  Power  Spectral  Density  of  the  Radiated  Noise 

After  the  substitution  of  Eq.  (15)  into  Eq.  (12)  we  find  that  the  radiated 
acoustic  pressure  is  of  the  form: 


Ax 


P(r,t)  = 


cos9 

2nrc 


iCr^.nj.t-r/c)  a(xQ+nL)  dr^dr^ 


(17) 


where  the  origin  of  our  coordinate  system  has  been  displaced  x^-units  downstream  in 
order  to  simplify  the  variables  in  our  integrals.  The  one-sided  physical  spectrum 
can  be  calculated  from  the  Fourier  integral  [23]: 


G(r,f)  = 2 


<p(r,t)  p(r,t+i)>  eia)T  dl 


(18) 


where  <p(r,t)  p(r,t+T)>  is  the  autocorrelation  function  of  p(r,t),  f is  the  fre- 
quency, and  w = 2irf.  In  general. 


<gg'>  = <g(t)  g ( t+T ) > 


T 

* 

lim  ^ 

T-x» 

0 


g(t)  g ( t+T ) d t 


(19) 


Therefore,  for  p(r,t): 
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2 

cos  0 


<pp  > = “TTI 

4ir  r c 


cKr^.r^.t-r/c)  I(4»1,<(>3,t-r/c+T)> 


•oUq+T^)  dn1dn3  , 


(20) 


where  for  the  time  being,  we  will  not  show  the  limits  of  integration.  It  will 
prove  useful  to  change  the  <p ^variables  to  spatial-separation  variables 
K-  = ^ - r^,  such  that  Eq.  (20)  becomes: 


2« 
cos  0 


<PP  > = 222 

4tt  r c 


<i(n1,n3»t-r/c)  i(ri1+^1* n3+?3»t-r/c+x)> 


^(Xq+Hj^)  o(x()+n1+?1) dfi1dn3  d?1d?3 


(21) 


If  we  assume  I(xi>t)  to  be  stationary  in  time  and  homogeneous  in  x3>  Eq.  (21)  may 
be  written  without  any  loss  of  generality  as: 


<PP'> 


2Q 
cos  0 


.222 
4tt  r c 


<101^0,0)  I(n1+51,53.T)> 


^(xq+Hj^)  a(x()+n1+C1)  dr^  dn3  dC1dC3 


(22) 


(We  note  that  it  makes  no  difference,  in  a stationary  process,  when  one  begins 

• • 

analyzing  an  event,  whether  at  t = t or  at  t'  = t +r/c.).  Because  <II'>  = 


9 9 

-3  <II'>/3x  , the  power  spectrum  of  a process  composed  of  Delta  functions  is 


equivalent  to  u times  the  power  spectrum  of  the  same  process  composed  of  unit 


impulse  functions.  Thus,  Eq.  (18)  becomes 
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G(r,f) 


2 2 
u>  cos  8 


. 2 2 2 , 
2tt  r c J 


<i(nL,o,0)  i(n1+C1,C3,x)> 


iU)T 


•oCxq+Hj^)  0(xo+n1+C1)  e dn1dn3  d^d^  dx 


(23) 


It  is  now  assumed  that  <II'>  can  be  treated  in  much  the  same  way  as  the 
analogous  space-time  correlation  function  that  occurs  in  the  case  of  a homogeneous, 
fully-developed  turbulent  boundary  layer.  In  particular,  we  assume  that  <II’> 
can  be  separated  [24,  25],  i.e. , let 


<II’>  = R2(ni’?3) 


(24) 


The  lateral  term,  R^,  can  be  integrated  in  the  form: 

L„ 


R3(ni>?3)  d?3 


d£3  = 2L3(n3) 


(25a) 


-L_ 


where  L3is  a transition-width  integral  scale.  Chen  and  Thyson  [20]  discuss  a 


spanwise  length  scale  in  the  region  of  transition  where  the  bursts  begin  to 
overlap;  that  is,  where  Ax.  They  write  this  scale  in  the  form 


A = 1.5  Ax  tan  a.  Let  us  set  L3(Ax)  = A and  L3(0)  = 0.  It  would  then  seem  that 


L3^1^  = 1,5ni  tan  a 


- n-j/4  for  a = 9.6" 


(25b) 


j 
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It  can  be  shown,  through  a Taylor  series  expansion,  that  0 = Tt~tl  is  only 
weakly-dependent  on  its  argument  through  the  range  of  integration  required  in 
Eq.  (23).  In  particular. 


c(xQ+e)  = Axq  (1  - 5x^)  - Bxq  (1  - 2x^) 


where 


A = 0. 0288pu  Re 
0 u 


Re  = u-/v,  which  is  the  unit  Reynolds  number, 
u 0 


The  maximum  value  which  £ assumes  is  Ax,  the  streamwise  extent  of  the  transition 
zone.  For  incompressible  boundary-layer  transition.  Ax  can  be  estimated  from  the 
following  emperical  relationship  [20,  22]  between  the  Reynolds  number  based  on  the 
transition  point,  x^,  and  that  based  on  Ax: 


Re.  = 60Re 
Ax  x. 
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In  connection  with  the  discussion  following  Eq.  (6),  we  would  expect  the  first  term 
(the  turbulent  wall  shear  stress)  of  Eq.  (26)  to  dominate,  so  that 

-1/5  -1/3 

O (x+e  ) = Ax  (l-12Re  ) 

0 max  0 Xq 

= o(x0)  . (28) 

This  approximation  implies  that  if  the  boundary  layer  were  completely  laminar  or 
completely  turbulent  over  the  distance  Ax,  then  the  mean  value  of  the  wall  shear 
stress  changes  insignificantly  over  this  streamwise  distance.  It  is  analogous  to 
the  parallel-flow  assumption  commonly  used  in  analyses  related  to  fully-developed 
turbulent  boundary- layer  flow. 

Differentiating  Eq.  (23)  with  respect  to  x^,  and  making  use  of  Eq's.  (25)  and 
(28),  yields: 


3G(r,f)  s 
3x_ 


(a  2cos26o2  (xq) 

.222 

4ir  r c 


niRi(,v^i’T> eiWT  dnl  d^i  dT 


(29) 


which  is  the  power  spectral  density  of  the  sound  radiated  per  unit  spanwise  width 
of  boundary-layer  transition. 

1.  Longitudinal  Space-Time  Correlation  Function 

The  longitudinal  space-time  correlation  function,  R^(n^,C^,T),  has  not  been 
previously  investigated  in  enough  detail  for  us  to  write  down  an  explicit  expression 
for  it;  however,  Shivitz  [26]  presented  a very  limited  set  of  data  for  R^  at 

= 0.6Ax.  Although  somewhat  heuristic  we  will  attempt  to  formulate  a model  for 
R^  based  on  our  physical  notions  of  the  transition  process.  In  regards  to  the 
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T-dependence,  we  let  the  random  sequence  of  Impulse  functions  that  make  up  our 
indicator  function  be  Poisson  distributed.  We  select  this  distribution  because 
it  is  a discrete  distribution  (which  it  necessarily  must  be)  for  a number  of  events 
(turbulent  bursts)  that  all  happen  at  random  times  with  an  average  of  N events  per 
unit  time  (our  burst  frequency).  Indeed,  the  experimental  investigation  of 
turbulence  bursts  by  DeMetz,  et.al.  [27]  does  suggest  a Poisson  process.  Under 
this  assumption,  the  classical  correlation  function  for  a random  telegraph  signal 
may  be  used  (see  Rice  [28]),  which  is  of  the  form  exp(-2N|r|). 

The  non-homogeneity  of  in  r)^  necessitates  weighing  it  by  the  standard 
deviation: 

<I2(n1,0)>  = Y(ni)  , (30) 

which  is  the  intermittency  factor  [from  Eq.  (13)].  At  a given  location,  r|^»  we 

expect  that  the  correlation  should  peak  when  T ■ £]/uc  represents  the  lag 

time  for  a given  turbulence  burst  to  convect  over  the  distance  at  a convection 

speed  u . However,  in  a frame  of  reference  moving  at  the  convection  speed  we 
c 

would  experience  a so-called  "moving-axis  decorrelation  [29]"  which  accounts  for 
the  fact  that  a burst  cannot  remain  perfectly  correlated  as  the  distance  increases 
beyond  some  characteristic  scale.  With  these  notions,  Eq.  (31)  appears  to  satisfy 
what  is  required: 

-0*15,1  -2N|t-£  /u 

R1(n1,C1,T)  * Y(nx)  e 1 e 1 (31) 

* 

where  a and  N are  presumed  to  depend  on  f|^* 
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By  way  of  a simple  experiment,  a pinhole  microphone  assembly  was  flush  mounted 
in  the  wall  of  a subsonic  wind  tunnel.  The  tunnel  speed  was  adjusted  so  that  the 
probe  was  in  the  transition  zone  of  the  wall  boundary  layer.  An  intermittency 
detector  enabled  us  to  generate,  electronically,  the  indicator  function  and  to 
determine  y and  N.  A real-time  correlation  analyzer  was  used  to  compute  the  auto- 
correlation function  of  I(x^,t).  The  result  of  this  experiment  is  shown  in  Fig.  6 
along  with  the  prediction  given  by  Eq.  (31)  for  = 0.  The  agreement  is  seen  to 
be  quite  good. 

With  regard  to  the  y and  N distributions,  we  note  that  Emmons  [9]  first  discussed 
them  and  developed  a probabilistic  model  to  predict  them.  He  assumed  the  existence 
of  a source-rate  density  function,  g(xQ,Zg,t),  which  specifies  the  rate  of  production 
of  turbulent  point-source  bursts  per  unit  area  on  the  surface  at  position  x^  = x^, 
x^  = Zq,  and  time  t^.  In  his  example,  Emmons  assumed  gCx^.z^.t^)  to  be  constant; 
however,  a later  investigation  by  Narasimha  [30]  which  made  use  of  the  data  of 
Ref.  [10],  showed  that  gCx^.z^.t^)  is  more  accurately  described  by  a Delta  function. 

In  particular,  g(x^)  = n<S  (x^-x^) , where  n is  defined  as  the  number  of  sources  per 
unit  length  per  unit  time  along  the  line  x^  = x^.  Farabee,  et.al  [11]  made  use  of 
this  form  of  g in  deriving  emperical  relations  for  y and  N.  They  are: 


Y(Z)  = 1-e 


-4.185Z 


(32) 


and 


N(Z)  = 1.272  ^ Ze_4,185Z 


(33) 


where  Z = r^/Ax.  Figure  7 shows  plots  of  N and  y as  computed  from  Eq's.  (32)  and 
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The  flare  parameter,  a , in  Eq.  (31)  is  considerably  more  difficult  to  predict 

because  of  the  lack  of  experimental  data.  Under  a fully-developed  turbulent  boundary 

* 

layer,  we  can  estimate  a from  Fig.  7 of  Blake’s  [29]  wall  pressure  data.  The 

* * 

smooth-wall  decay  envelope  can  be  approximated  by  exp (-£^/26  ),  where  6 is  the 

* * -1 

displacement  thickness.  Thus,  we  let  a (Z 2 1)  = (26  ) . Within  the  transition 

zone,  we  can  make  use  of  the  space-time  correlations  measured  by  Shivitz,  ibid, 

* 

Fig.  24.  Here,  Z = 0.6  and  the  decay  envelope  can  be  approximated  by  exp (-^/HOA  ), 
* 

where  the  6 is  that  measured  shortly  downstream  of  where  the  flow  became  fully 
* * -1 

turbulent;  hence,  a (Z  = 0.6)  = (1106  ) . Farther  upstream  in  the  transition 

* 

zone,  we  can  only  speculate  as  to  what  value  a may  have.  Because  there  will  be 

fewer  bursts  of  turbulence  as  Z -*■  0,  we  might  anticipate  that  the  moving-axis 

* 

decorrelation  becomes  less  severe,  meaning  a should  decrease  as  Z -*•  0.  We  will 

* 

assume  this  to  be  the  case,  and  also,  that  a can  never  become  smaller  than 

-1  * -1  * 

(Ax)  ; thus,  we  let  a (Z  = 0)  = (Ax)  . Using  the  three  values  of  a discussed 

* 

above,  and  the  values  of  6 and  Ax  measured  by  Shivitz,  the  following  equation  for 
* 

a can  be  constructed: 


a*  = (1+83.35Z8)/Ax  , (34) 

which  is  shown  graphically  in  Fig.  8. 

The  space-time  correlation  function  given  by  Eq.  (31)  was  computed  for  Z = 0.6 
and  uq/Ax  = 467  Hz  (typical  value  for  water) . Figure  9 shows  the  resulting  three- 


dimensional surface. 
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L: 


2.  Integration  of  Eq.  (29) 


If  we  now  substitute  Eq.  (31)  into  Eq.  (29),  and  change  the  variable  (t-£,/u  ) 

1 c 


to  T ' , we  then  have : 


Ax  Ax-ru 


9G(r,f) 

3x3 


2 2 2 
to  cos  00 


.222 

4tt  r c 


'liY 


ICjl  “A 


-“b'l  im-  . , „r  . 

i e dt'  d^  dnx 


o -nn 


2 2 

where  = to/uc«  The  last  integral  is  well-known  and  equal  to  2N/ (4N  +to  ). 


Because  the  power  spectrum  must  be  real,  the  above  equation  reduces  to: 


Ax 


9G(r,f)  „ to2cos20q2 


3x. 


.222 

2tt  r c 


nx  Y N ] 


2 2 
4lT+w 


* 

a C, 


Ax-n 


e cos  k^d^  + 


1 * 


e cos  kc?1d^1  , dT^OS) 


rni 


The  integrations  over  can  be  performed.  Without  showing  the  details,  but  noting 


that  the  following  identity  was  used : 


-1  b. 


a sin0  + b cos0  = /a2+b2  sin(0+tan  — ) , 

a 


and  that  the  change  of  variable  from  ^ to  Z = P^/Ax  was  also  made,  Eq.  (35)  becomes: 


ac(r.f)  o2(xo)cos20uo(Ax)2 

-*£-  = TTT2 F(kcAx,a  ^/uQ)  * 

3 2tt  r c 


(36) 


. 
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with 


* . 2 
F(k  Ax, a Ax,u  / u_)  = (k  Ax) 
c cue 


,2  . -4.185Z2  -4.185Z2 


6.472 

V 

2 2 
,2  -8.37Z * . .2 

Z e +(k  Ax) 

/(a  Ax)2+(kcAx)2 

CJ 

2a  Ax  j - e'a  AxZsin(kcAxZ-<J))  + e_0t  Ax(1"Z)sin[kcAx(l-Z)-<|>]  ) dZ  , (37) 


/(a  ' Ax) 2 +(k  Ax) 
c 


—1  * 

and  = tan  (a  Ax/kcAx).  In  the  development  of  Eq.  (36),  it  has  been  assumed  that 
u c/ Uq  is  independent  of  frequency. 

Equation  (37)  cannot  be  integrated  conveniently;  however,  because  of  its  non- 

dimensional  form,  it  need  only  be  integrated  once  for  a few  discrete  values  of 

* 

uc'u0  an<*  a ran8e  k£Ax  [assuming  Eq.  (34)  describes  the  behavior  of  a Ax].  This 

integration  has  been  performed  for  U(;/Uq  = 0.6  and  0.8  by  Simpson's  rule.  The  result 

is  presented  graphically  in  Fig.  10.  We  see  that  boundary-layer  transition  generates 

significant  noise  only  in  a frequency  band  centered  about  f = u /2irAx. 

o c 

2 

Because  o(Xg)  is  proportional  to  pu^  , Eq.  (36)  suggests  that  the  spectral  level 

increases  according  to  for  constant  Ax.  However,  Eq.  (27)  suggests  that  Ax  is 

3 

inversely  proportional  to  Ug,  so  Eq.  (36)  predicts  a Uq  velocity  dependence  even 

though  the  source  of  noise  is  dipole.  The  reason  for  this  unexpected  behavior  is 

that  the  solution  is  derived  in  terms  of  the  detailed  hydrodynamics.  These  details 

result  in  a source  area  that  changes  with  velocity.'  Thus,  to  be  consistent  with 

less-rigorous,  flow  noise  order  of  magnitude  estimates,  the  radiation-field  pressure 

should  be  normalized  by  the  area  of  the  source  region,  which  is  proportional  to 

Ax.  The  mean-square  value  of  the  sound  pressure  radiated  per  unit  area  would  then 

2 

be  given  by  the  integration  of  Eq.  (36)  over  frequency  divided  by  (Ax)  . A 5th 
power  velocity  dependence  would  thus  be  retained. 
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E.  Radiation  Efficiency 

The  radiation  efficiency,  n > is  defined  as  the  ratio  of  acoustic  power,  N , 

tr  • a 

radiated  to  the  farf ield-to-the  work  expended  per  unit  time  by  the  hydrodynamic 
motions  within  the  acoustic  source  region,  Nh»  By  re-writing  Eq.  (36)  in  the  form: 


3G(r ,f ) 
3x„ 


2a  s 

cos  0 <Pq  > 


we  see  that 


2 

<pn  > 

II  II 

N = 

— y — du 

a 

4 

pc 

4 

2 2 

cos  0r  sin0d0d(J) 


(38) 


0 0 


Upon  substituting  for  <pQ  >,  Eq.  (38)  reduces  to: 


N 

a 


2 2 

2o  Uq  Ax 
3irpc^ 


(39) 


where 


00 

* 

F(k  Ax, a Ax,u  /u  ) d(k  Ax) 
c c u c 

4 

o 


The  magnitude  of  F is  denoted  in  Fig.  10  for  each  of  the  two  convection  speed  ratios 

considered.  Essentially,  F is  independent  of  u /u„  and  is  of  order  0.1. 

c 0 
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Because  Eq.  (39)  is  the  acoustic  power  radiated  per  unit  spanwise  width  of 
transitional  flow,  it  is  necessary  to  compute  the  hydrodynamic  power  also  on  a per 
unit  width  basis.  The  energy  expended  per  unit  time  per  unit  width  by  transitional 
flow  may  be  expressed  by: 


x0+Ax 


Nh  = u0 


To  dxi 


(40) 


where  tq  is  the  mean  value  of  the  wall  shear  stress.  Because  bursting  flow  occurs 
between  the  limits  of  integration,  we  would  expect  that 

• To  (xl>  = TT  <xiMxi>  + [l-y(x1)  ] tl  (x^  . (41) 

Because 


To  * (TT  - Tl>  Y + TL 


where  (tt  - T^)  = O ~ 0.9tt>  we  let 


Tq  (xx)  = y(.Xj)  o (xx) 


(42) 


The  intermittency  distribution  can  be  calculated  using  Eq.  (32).  Again  we  use 
Eq.  (28)  to  arrive  at 

* 

1 


Nh  = uqAxcXx0) 


1 - 


a"4-1852  dZ 
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Through  a change  of  integration  variable,  the  above  integral  can  be  expressed  in 
terms  of  the  error  function,  i.e.. 


N,  = u Axa[l-0.488  ^ erf (2. 048)] 

ho  2 


= 0.572  UqAx  a 


(43) 


I 

i 


The  radiation  efficiency  is  now  obtained  by  dividing  Eq.  (39)  by  Eq.  (43), 
which  results  in 


0.37 


OU0F 

pc3 


(44) 


Again  noting  that  o-  pu^  , we  see  that  the  radiation  efficiency  is  proportional  to 
the  cube  of  the  free-stream  Mach  number.  This  is  analogous  to  what  Landahl  [7] 
found  for  the  shear  stress  contribution  to  fully-developed  turbulent  boundary-layer 
radiated  noise,  i.e.. 


(45) 


where  u^  is  the  friction  velocity  and  = u^/c. 

We  would  like  to  compare  the  radiation  efficiency  of  transitional  flow  to  that 
of  a fully-developed  turbulent  boundary  layer  (TBL)  flow.  Equation  (45)  cannot  be 
used  in  this  comparison  due  to  the  lack  of  a proportionality  factor.  An  analysis 
given  by  Tam  [31],  however,  may  be  used  in  this  regard.  Tam  did  not  derive  an 
explicit  expression  for  the  efficiency,  but  did  give  an  expression  for  the  acoustic 
power  generated  per  unit  surface  area  of  boundary  layer  turbulence.  In  particular. 


r 
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4itt„ 


N. 


TBL  pc 


F(S,M)  dS 


(A6) 


where  the  integral  is  presented  graphically  in  Reference  [31]  as  a function  of  Mach 
number,  M.  Tam  concluded  that  the  magnitude  of  this  integral  increases  rather 
rapidly  with  Mach  number  (slightly  faster  than  M2).  We  have  thus  been  able  to 
approximate  it  by: 


F(S,M)  dS  - 0.015M2  (M  > 0) 


(47) 


The  hydrodynamic  power  generated  per  unit  surface  area  of  boundary  layer  turbulence 
is  given  by  u^x^,  so  we  would  expect  that: 


'TBL 


0.188  XT  M2 
Pcun 


(48) 


By  letting  a - X , we  find  that 


\r.  ~ 1,97  F nTBL  ~ 0,2  nTBL 


(49) 


This  result  implies  that  the  wall  shear  stress  fluctuations  that  occur  in  inter- 
mittent transition  zone  flow  generate  noise  less  efficiently  than  do  the  super- 
sonic Fourier  components  of  the  normal  stress  fluctuations  [31]  in  a fully— developed 


turbulent  boundary  layer. 
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III.  EXAMPLE 

In  order  to  verify  the  analysis  presented  In  this  paper,  one  must  conduct  a very 
sophisticated  experiment.  This  author  envisions  such  an  experiment  as  being  composed 
of  a large,  nearly  rigid,  thin  flat  plate  mounted  in  an  anecholc  wind  tunnel  (air 
is  better  than  water  because  of  the  rigidity  requirement) . The  span  of  the  plate 
must  be  very  large  as  to  minimize  the  edge  noise  received,  say  at  the  center  of  the 
plate.  Because  the  theory  predicts  a cosine  directivity  curve,  one  must  use  flush 
mounted  microphones  in  order  to  receive  the  peak  sound  level.  They  should  be  placed 
in  the  laminar  flow  regime  of  the  plate's  boundary  layer  as  to  minimize  near-field 
pressure  fluctuation  effects.  Perhaps  most  importantly,  the  operating  Reynolds 
number  should  be  adjusted  so  that  fully-developed  turbulent  flow  occurs  very  near 
the  trailing  edge  of  the  plate  (x^  + Ax  = length  of  plate) . The  trailing  edge  must 
be  designed  to  minimize,  as  much  as  possible,  any  unsteady  wake  flow.  Unfortunately, 
an  experiment  such  as  this  one  has  not  yet  been  carried  out. 

Nevertheless,  it  is  tempting  at  this  point,  to  make  a prediction  using  the 
subject  theory  for  a situation  in  which  noise  was  measured  flush  to  the  surface  and 
within  the  laminar  boundary  layer  of  a test  vehicle.  We  choose  those  data  presented 
by  Nisewanger  and  Sperling  [32].  In  their  experiment,  a buoyantly-propelled  axisym- 
metric  body  was  instrumented  with  flush-mounted  hydrophones  and  noise  spectra  at 
various  locations  on  the  surface  were  measured  in  fresh  water  at  a vehicle  speed  of 
19.4  m/s.  The  transition  point  was  found  to  be  at  an  arc  length  distance  of  16  cm. 
Equation  (27)  is  used  to  calculate  Ax;  we  find  it  to  be  5.8  cm  using  the  kinematic 
viscosity  for  fresh  water  at  21°C.  Measurement  station  no.  4 of  Ref.  32  is  selected 
because  it  is  in  the  laminar  flow  region,  2.6  cm  forward  of  the  transition  point. 
Letting  9=0°,  u^/uq  ■ 0.8,  and  multiplying  Eq.  (36)  by  the  circumference  of  the 
test  body,  the  resulting  comparison  of  theory  and  experiment  is  shown  in  Fig.  11. 
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As  might  be  expected,  the  agreement  is  not  very  good.  The  differences  in  level 
may  be  attributable  to  several  reasons.  Because  the  test  vehicle  supports  a substantial 
area  of  fully-developed  turbulent  boundary  layer  flow,  the  data  may  be  contaminated 
by  this  source  of  noise.  Also,  as  noted  by  Nisewanger  and  Sperling,  approximately 
50  percent  of  the  measured  flow  noise  is  due  to  shell  vibrations.  From  the  more 
fundamental  point  of  view,  the  shear  stress  mechanism  of  sound  generation  by  a 
transition  zone  may  not  be  the  only  mechanism  involved.  We  neglected  the  Reynolds 
stress  contributions  through  an  argument  of  relative  efficiency  between  dipole  and 
quadrupole  sources.  It  would  seem  reasonable  that  the  fluctuating  shear  stress, 
being  a viscous  quantity,  may  be  considerably  smaller  than  the  fluctuating  Reynolds 
stress;  perhaps  even  small  enough,  that  the  relative  difference  in  radiation 
efficiency  is  out-weighed. 

IV.  CONCLUSIONS 

We  have  presented  an  analysis  of  the  radiated  sound  due  to  boundary-layer 
transition.  Our  principle  assumptions  were  that  the  surface  is  infinite  in  extent, 
planar,  and  rigid,  that  the  transition  process  includes  a finite  region  intermittent 
flow  where  the  boundary  layer  fluctuates  randomly  between  laminar  and  turbulent,  and 
that  it  is  this  intermittent  flow  that  generates  the  noise.  Through  use  of  Lighthill's 
analogy,  we  showed  that  the  fluctuations  between  laminar  and  turbulent  boundary- 
layer  flow  give  rise  to  dipole,  quardupole,  and  octupole  noise  sources.  On  the 
basis  of  a very  low  Mach  number  assumption,  we  treated  only  the  dipole  contribution 
in  detail;  this  contribution  being  due  to  the  fluctuating  wall  shear  stresses.  We 
assumed  that  the  flow  is  statistically  homogeneous  in  the  spanwise  directions, 
non-homogeneous  in  the  streamwise  direction,  and  stationary,  but  Poisson  distributed 
in  time.  The  power  spectral  density  for  the  acoustic  pressure  radiated  per  unit 
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spanwise  width  was  derived.  The  level  of  the  mean-square  radiated  pressure  was  found 
to  depend  upon  the  square  of  the  difference  between  the  turbulent  wall  shear  stress 
and  the  laminar  wall  shear  stress  that  would  occur  at  the  beginning  of  transition, 
upon  the  square  of  the  streamwise  distance  over  which  turbulent  bursts  occur,  and 
upon  the  free-stream  velocity.  The  spectrum  rises  at  6 dB/octave  at  low  frequencies 
and  drops  off  at  6 dB/octave  at  high  frequencies.  A broad  peak  is  centered  about 
kcAx  = 1.0.  An  expression  for  the  radiation  efficiency  was  derived  and  compared  to 
the  radiation  efficiency  of  a fully-developed  turbulent  boundary  layer.  This 
comparison  showed  that  transition  generates  noise  at  approximately  20  percent  the 
efficiency  of  fully-developed  turbulent  boundary  layer  flow. 

We  compared  a prediction  using  this  theory  with  the  noise  spectrum  measured 
near  the  laminar-to-turbulent  transition  region  of  a buoyantly-propelled  flow  noise 
research  vehicle.  The  agreement  was  considered  to  be  poor.  Based  on  this 
comparison  and  on  the  magnitude  of  the  radiation  efficiency  it  can  be  concluded  that 
the  shear  stress  fluctuations  that  occur  in  natural  boundary-layer  transition  are 
not  strong  sound  radiators.  However,  boundary-layer  transition  is  a very  complicated 
hydrodynamic  phenomenon,  and  the  neglect  of  fluctuating  Reynolds  stresses  on  the 
supposition  that  they  give  rise  to  inefficient  quardupoles  and  octupoles  may  not  be 
fully  justified.  A detailed  examination  of  the  first  term  of  Eq.  (9)  appears  to  be 
required  but  the  modeling  of  u^u^  i*1  an  intermittent  flow  regime  may  prove  to  be 
very  difficult.  In  order  to  simplify  this  task,  it  may  prove  expedient  to  first 
investigate  the  third  and  sixth  terms  of  Eq.  (10)  which  describe  the  normal  stress 
contributions.  Such  analysis  must  be  performed  in  wavenumber/frequency  space  so 
that  those  Fourier  components  with  phase  velocities  (u>/k)  greater  than  or  equal  to 
the  velocity  of  sound  can  be  identified  (if  present).  It  would  be  these  components 
that  lead  to  sound  radiation;  any  others  would  cause  near-field  pressure  fluctuations. 
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disprove  the  notion  that  boundary-layer  transition  is  an  important  source  of  low 
Mach  number  flow  noise.  The  space-time  correlation  functions  at  various  locations 
within  a transition  zone  should  be  measured.  These  functions  would  not  only  aid 
in  the  modeling  presented  in  this  paper s but  would  provide  the  groundwork  necessary 
to  proceed  with  further  analysis  such  as  that  suggested  above. 
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Figure  3.  (a)  A schematic  representation  of 

the  image  flow  concept;  and, 

(b)  the  way  in  which  multipole 
sources  in  the  real  flow  combine 
with  e<,  '.valent  sources  in  the 
image  flow. 
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Figure  4.  The  indicator  function  composed  of  unit  impulse 
functions  and  its  first  time  derivative. 
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Figure  5.  A schematic  representation  of  I(x^,t)  in  space  and 
time,  where  y is  its  time  average  value  and  N is 
the  expected  number  of  bursts  per  unit  time. 
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Figure  11.  Comparison  of  the  theory  with  experimental  spectra  measured  on  the 
surface  of  a buoyantly-propelled  flow  noise  research  vehicle  [32]. 
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